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We study finite size effects in superconducting metallic grains and determine the BCS order 
parameter and the low energy excitation spectrum in terms of size, and shape of the grain. Our 
approach combines the BCS self-consistency condition, a semiclassical expansion for the spectral 
density and interaction matrix elements, and corrections to the BCS mean-field. In chaotic grains 
mesoscopic fluctuations of the matrix elements lead to a smooth dependence of the order parameter 
on the excitation energy. In the integrable case we observe shell effects when e.g. a small change in 
the electron number leads to large changes in the energy gap. 
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PACS numbers: 74.20.Fg, 75.10.Jm, 71.10.Li, 73. 21. La 

Since experiments by Ralph, Black, and Tinkham [l[ 
on Al nanograins in mid nineties, there has been con- 
siderable interest in the theory of ultrasmall supercon- 
ductors (see @, 13] for earlier studies). In particular, 
finite-size corrections to the predictions of the Bardeen, 
Cooper, and Schriffer (BCS) theory for bulk supercon- 
ductors [H have been studied joT-flO] within the exactly 
solvable Richardson model [ll|. Pairing in specific po- 
tentials, such as a harmonic oscillator potential 12j and 
a rectangular box, 13, 141 and mesoscopic fluctuations 
of the energy gap [la . [l6[ have been explored as well. 
Nevertheless, a comprehensive theoretical description of 
the combined effect of discrete energy spectrum and fluc- 
tuating interaction matrix elements has not yet emerged. 
We note that the Richardson model alone cannot provide 
such a description as it does not allow for mesoscopic fluc- 
tuations of the matrix elements. 



In the present paper we develop a framework based on 
the BCS theory and semiclassical techniques that permits 
a systematic analytical evaluation of the low energy spec- 
tral properties of superconducting nanograins in terms of 
their size and shape. Leading finite size corrections to 
the BCS mean-field can also be taken into account in our 
approach. Our main results are as follows. For chaotic 
grains, we show that the order parameter is energy de- 
pendent. The energy dependence is universal, i.e. its 
functional form is the same for all chaotic grains. The 
matrix elements are responsible for most of the deviation 
from the bulk limit. In integrable grains, we find that the 
superconducting gap is strongly sensitive to shell effects, 
namely, a small modification of the grain size or number 
of electrons can substantially affect its value. 



We start with the BCS Hamiltonian, H = 
E„ CT CnC&trCn* ~ E„,„' T n,n' cLcLc n ' ;<v T , where c na an- 



nihilates an electron of spin a in state n, 



I(e n ,e n ,)=\V6 / ^l(r)vbl,(r)dV 



(1) 



are matrix elements of a short-range electron-electron in- 
teraction, A is the BCS coupling constant, and vb n and 
e n are eigenstates and eigenvalues of the one-body mean- 
field Hamiltonian of a free particle of mass m in a clean 
grain of volume V. Eigenvalues e n are measured from the 
Fermi level ep and the mean level spacing S = l/iA r p(0), 

where ^tf(O) = 2^j- (^f) 3 2 \ff~F is the spectral density 
at the Fermi level in the Thomas-Fermi approximation. 

Our general strategy can be summarized as follows: a) 
use semiclassical techniques to compute the spectral den- 
sity v{e) — J2 n — e «) ano - -^( e ' e ') as ser i es i n a small 
parameter 1/kpL, where Uf is the Fermi wavevector and 
L ~ I/ 1 / 3 is the size of the grain b) solve the BCS gap 
equation in orders in 1/kpL c) evaluate the low energy 
spectral properties of the grain such as the energy gap, 
excitation energies, and Matveev-Larkin parameter [fj( 
including finite size corrections to the BCS mean-field. 
The results thus obtained are strictly valid in the region, 
kpL 1 (limit of validity of the semiclassical approxi- 
mation), S/Aq < 1 (limit of validity of the BCS theory), 
and I 3> £ S> L (condition of quantum coherence) where 
£ = Hvf/ A-o is the superconducting coherence length, vf 
is the Fermi velocity, / is the coherence length of the sin- 
gle particle problem and A is the bulk gap. We note 
that in Al grains [§] £ w 1600nm and I > lOOOOnm for 
temperatures T < AK. Therefore the region I ^> £ 3> L 
is accessible to experiments. 

Since the matrix elements I(e, e') are energy dependent 
the BCS order parameter A(e) also depends on energy. 
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The self-consistency equation for A(e) reads 
" E ° A( e ')J(e, e ') 



A(e) 



v(e')de', (2) 



-E D 2 Je' 2 + A(e') 2 



where Ep, is the Debye energy. In the limit V — > oo, 
the spectral density in the 2i?£> energy window near 
the Fermi level can be taken to be energy indepen- 
dent and given by the Thomas-Fermi approximation, 
v(e) — zatf(O), matrix elements are also energy indepen- 
dent, I(e, e') = XS, and the gap is equal to its bulk value, 
Ao = 2Ep,cT x . As the volume of the grain decreases the 
mean level spacing increases and eventually both v(e) 
and /(e, e') deviate from the bulk limit. 

Semiclassical evaluation ofu(e). The spectral density 
in a 3d grain, 



p(e')~v TF (0) [1 + 5(0) 



(3) 
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consists of a monotonous part, g(0) — _i_ 4fcf y 
and an oscillatory contribution gi(e'). Here S and 
C denote the surface area and mean curvature of the 
grain, respectively, and upper/lower signs stand for Neu- 
mann/Dirichlet boundary conditions. The oscillatory 
contribution, to leading order, is given by the Gutzwiller 
trace formula 13, HI, 



I 



klv^ Ape 



e'k F L p 
2e F 



(4) 



where both the amplitude A p and the topological in- 
dex (3 P depend on classical quantities only [18j| . The 
summation is over a set of classical periodic orbits p of 
length L p . For isolated grains Dirichlct is the most nat- 
ural choice, but we also include Neumann to illustrate 
the dependence of our results on boundary conditions. 
Only orbits shorter than the quantum coherence length 
I of the single-particle problem are included. This effec- 
tively accounts for inelastic scattering and other factors 
that destroy quantum coherence. Here we focus on the 
limit / 3> £, the case I ~ £ will be discussed elsewhere 
[20I ] . In Eq. (HJ) classical actions hk(e')L p are expanded 

as k(e') m kp + e'k F /2e F . The amplitude A p increases 

1/2 

by a factor (k F L) 1 3> 1 for each of the symmetry axes 
of the grain. 

Semiclassical evaluation of I(e,e'). For integrable sys- 
tems J(e, e') depends on details of the system. In a rect- 
angular box it is simply /(e, e') = X5 but in most other ge- 
ometries an explicit expression in terms of classical quan- 
tities is not available. In the chaotic case the situation is 
different. As a result of the quantum ergodicity theorem 
[1^ ^ i s wen justified to assume that for systems with 
time reversal symmetry (the only ones addressed in this 
paper), ^(r) = y(l + 0(l/k F L)). In order to explicitly 
determine deviations from the bulk limit we replace ip^ (r) 



in I n:n > with (ip 2 (r)) eri , where (. . .) c stands for an energy 
average around e. The single-particle probability density 
is thus effectively averaged over a small energy window 
resembling the effect of a finite coherence length. 
Substituting {ip 2 (f)) e into I n ,n', we obtain 
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+ I(e F ,e,e') 



where 



I(e F , e, e') = / short (e F ) + I lo ^(e F , e - e') 



(5) 



(6) 



can be split into two parts coming from short and long 
orbits. Short orbits involve a single reflection at the grain 
boundary and result in a monotonous contribution 



short 



irS 
AkpV" 



(7) 



while the contribution of long orbits depends on the en- 
ergy difference e — e' 



I i ^(e F ,e-e') = ^U l 



(8) 



with Hi(w) = J ^ ,< D 2 cos [wkpL^dr, where the sum 
is over all non-zero classical paths (not periodic orbits) 
7(7") starting and ending at a given point r inside the 
grain [2~0 ] and the amplitude D 1 is defined in Refs. [It], 
2oL 21 1 . The integral stands for an average over all points 
r inside the grain. The explicit evaluation of Tli(w) for 
a given geometry requires in principle the knowledge of 
all classical paths L 1 up to length I, However, for / ^> L, 
one can use a sum rule for classical closed orbits (23j to 
obtain 



sin(wk F I) 
wk F 



(9) 



Solution of the gap equation. First, let us consider 
chaotic grains. Here we present only the final answer for 
the 3d case deferring a more detailed account, including 
the 2d case, to Ref. poj . Writing the gap function A(e) 
formally as a series in l/k F L, 



A(e) = A fl + /W+/( 2 )+/( 3 )( £ ) 



(10) 



substituting it into Eq. ([2]), and using the above expres- 
sions for the density of states and interaction matrix el- 
ements, we derive 



/ 



(i) _ 



1±1 ttS 

A Ak F V' 



(11) 



where ± stands for Neumann (+) and Dirichlet (— ) 
boundary conditions. Note that to leading order the 
combined effect of the interaction matrix elements and 
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FIG. 1: Superconducting order parameter A(e) in units of 
the bulk gap Ao for chaotic Al grains (/cf = 17.5nm _1 ,<5 = 
7279/iV, Ao ~ 0.24meV) as a function of energy e counted 
from the Fermi level. Different curves correspond to grain 
sizes (top to bottom) L = 6nm,kpL = 105,5/Ao = 
0.77) (Dirichlet and Neumann boundary conditions), L — 
8nm,k F L = 140,<5/A = 0.32 (Dirichlet), and L = 
Wnm,k F L = 175, S/A = 0.08 (Dirichlet). The leading con- 
tribution comes from the energy dependent matrix elements 
I(e,e') given by Eq. 



the density of states have very different consequences on 
the gap, depending on the kind of boundary conditions. 
For Dirichlet the leading finite size corrections to the gap 
vanishes. 

The second order (l/(kpL) 2 ) correction reads 



Xf {2) 
where, 



2C 



■5(0), 



that mesoscopic corrections given by Eqs. (|11I12[) and 
(fl"4)) always enhance A(0) as compared to the bulk value 
Ao- Fig. Q] shows the gap function A(e) for Al grains of 
different sizes L, where we used (see [1() kp « 17.5nm _1 , 
A w 0.18, and 5 w 7279/N meV, where N the number of 
particles. 

Several remarks are in order: a) the smoothing of the 
spectral density energy dependence in Eq. (fT2|) caused 
by a cutoff function W is a superconductivity effect not 
related to the destruction of quantum coherence, b) the 
energy dependence of the gap is universal in the sense 
that it does not depend on specific grain details, c) the 
matrix elements /(e, e') play a crucial role, e.g. they are 
responsible for most of the deviation from the bulk limit 
in Fig. [TJ d) the requirement £ 3> L used to derive Eq. (|9|) 
is well justified for nanograins since L ~ lOnm, while 
£ - 10 4 nm. 

We now turn to the integrable case. Probably the 
simplest example is that of a rectangular box, since 
in this case the interaction matrix elements are simply 
I(e, e') = A<5. The calculation is simplified as now the 
order parameter is energy independent. We have 



A = Ao [1 + + fW) + f 

where f( n > cx (fcp£) -n A -1 . We obtain 
A/ (1) =5(0)+3 (1) (0), 



f(2) 



(15) 



A/ 



(3/2) 



E 



-(3/2) 



(0), 



(16) 



A/ (2) = £3^(0) + 2(0)], 



0(0) = 



2vr 



A p W{L p /0 cos(fc F L p + f3 p ) (12) 



and, 



W(L p /0 = $ 



dt 



cos(L p t/£) 



(13) 



exponentially suppresses periodic orbits longer than £. 

The third order correction (included in the definition 
of S) is energy dependent, 



/ (3) (e) 



7rA<5 



A,, 



A 2 
^0 



(14) 



Note that a) 8/ Ao <SC 1 is an additional expansion param- 
eter, therefore the contribution fT4"|) can be comparable 
to lower orders in the expansion in 1/kpL and b) the or- 
der parameter A(e) has a maximum at the Fermi energy 
(e = 0) and slowly decreases on an energy scale e ~ Ao as 
one moves away from the Fermi level. One can also show 



where cx (kpL)~ k denotes the oscillating part of the 
spectral density and indexes i and j take values 1, 2, 
and 3 in three dimensions. Explicit expressions for g^ k \ 
g!- k , and g^ in terms of periodic orbits for a rectangu- 
lar box can be found in Ref. [l8| (the cutoff function in 
our case is given by Eq. (jT3j)) . We note that: a) Eq. (fT5l) 
is also obtained by expanding the standard expression 
of the bulk gap A = 2Ep> exp(— j/tf(0)/^(0)A) in pow- 
ers of (kpLX) -1 with v(0) given by Eq.Q. b) unlike 
the chaotic case, the leading smooth correction to the 
bulk limit does not vanish for any boundary condition, 
c) smooth and oscillating corrections are of comparable 
magnitudes. 

Shell effects and fluctuations. Motivated by previous 
studies for other fermionic systems such as nuclei and 
atomic clusters (see e.g. Ref. H^), we investigate shell 
effects in metallic nanograins. In particular, we are inter- 
ested in the fluctuations of the BCS gap with the num- 
ber of electrons on the grain. As an illustration let us 
consider a cubic geometry. To determine the gap, we 
solve the gap equation ^ numerically and determine the 
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FIG. 2: Superconducting order parameter A in units of the 
bulk gap Ao for a cubic Al grain as a function of the ratio 
Ao/5, where S is the mean level spacing. Black crosses cor- 
respond to the exact numerical solution of the gap equation 
([2]) , while the red circles represent the semiclassical analytical 
expression (JT5J ■ 



Fermi energy for a given number of electrons TV by in- 
verting the relation 2 J F v(e)de — N . We find a good 
agreement between numerical results and the semiclassi- 
cal expansion (fl"5"|) , see Fig. [5J We also observe that a 
slight modification of the grain size (or equivalently the 
number of electrons N or the mean level spacing S) can 
result in substantial changes in the value of the gap, see 
Fig. O The typical magnitude of fluctuations of the gap, 

3^ ~ J [IH is consistent with our results (see Fig. 
2)°. 

Low energy excitations. Having solved the gap equa- 
tion (HI), one can evaluate low energy properties of the 
grain taking into account finite size corrections to the 
BCS mean-field approximation. For example, the energy 
cost for breaking a Cooper pair in an isolated grain is 

H, 



AE = 2A(0) - 5, 



(17) 



where A(0) is the solution of equation |(3J) taken at the 
Fermi energy and is given by Eqs. (fT0|) and (fTS"|) for 
chaotic and rectangular shapes, respectively. We note 
that the correction to the mean-field (—&) has been evalu- 
ated [13] for constant interaction matrix elements. Never- 
theless, since the deviation of matrix elements from a con- 
stant energy independent value is itself of order (fc^i)" 1 , 
Eq. (JTT|) is accurate up to terms of order (5/ A )(kpL)~ 1 , 
which are negligible as compared to the ones we kept in 
Eqs. CLZ]), (HOJ), and (fT5)l. 

Similarly, the Matveev-Larkin parity parameter [B[ 
reads A p = E 2N+1 - \{E 2N + E 2N+2 ) = A(0) - f, 
where En is the ground state energy for a superconduct- 
ing grain with N electrons. Quasiparticle energies are 

e 2 + A(e) 2 plus corrections to mean-field, which can 
be determined using the approach of Ref . [lfj ■ 

We see that finite size corrections to the BCS mean- 



field approximation are comparable to the energy depen- 
dent correction (|10[) obtained within mean-field, but have 
an opposite sign. We also note that our approach of ex- 
panding around the bulk BCS ground state is applicable 
only when 6 -C Ao, i.e. when corrections to the BCS 



mean- field approximation are small 25 1. 

To conclude, we have determined the low energy ex- 
citation spectrum for small superconducting grains as a 
function of their size and shape by combining the BCS 
mean-field, semiclassical techniques and leading correc- 
tions to the mean-field. For chaotic grains the non-trivial 
energy dependence of the interaction matrix elements 
leads to a universal smooth dependence fTJ| of the gap 
function on excitation energy. In the integrable case we 
found that small changes in the number of electrons can 
substantially modify the superconducting gap. 
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